Microbial induced wettability alteration with implications for Underground Hydrogen Storage.

Characterization of the microbial activity impacts on transport and storage of hydrogen is a crucial aspect of successful Underground Hydrogen Storage (UHS). Microbes can use hydrogen for their metabolism, which can then lead to formation of biofilms. Biofilms can potentially alter the wettability of the system and, consequently, impact the flow dynamics and trapping mechanisms in the reservoir. In this study, we investigate the impact of microbial activity on wettability of the hydrogen/brine/rock system, using the captive-bubble cell experimental approach. Apparent contact angles are measured for bubbles of pure hydrogen in contact with a solid surface inside a cell filled with living brine which contains sulphate reducing microbes. To investigate the impact of surface roughness, two different solid samples are used: a "rough" Bentheimer Sandstone sample and a "smooth" pure Quartz sample. It is found that, in systems where buoyancy and interfacial forces are the main acting forces, the impact of biofilm formation on the apparent contact angle highly depends on the surface roughness. For the "rough" Bentheimer sandstone, the apparent contact angle was unchanged by biofilm formation, while for the smooth pure Quartz sample the apparent contact angle decreased significantly, making the system more water-wet. This decrease in apparent contact angle is in contrast with an earlier study present in the literature where a significant increase in contact angle due to microbial activity was reported. The wettability of the biofilm is mainly determined by the consistency of the Extracellular Polymeric Substances (EPS) which depends on the growth conditions in the system. Therefore, to determine the impact of microbial activity on the wettability during UHS will require accurate replication of the reservoir conditions including surface roughness, chemical composition of the brine, the microbial community, as well as temperature, pressure and pH-value conditions.

Calculating Thermodynamic Factors for Diffusion Using the Continuous Fractional Component Monte Carlo Method.

Thermodynamic factors for diffusion connect the Fick and Maxwell-Stefan diffusion coefficients used to quantify mass transfer. Activity coefficient models or equations of state can be fitted to experimental or simulation data, from which thermodynamic factors can be obtained by differentiation. The accuracy of thermodynamic factors determined using indirect routes is dictated by the specific choice of an activity coefficient model or an equation of state. The Permuted Widom's Test Particle Insertion (PWTPI) method developed by Balaji et al. enables direct determination of thermodynamic factors in binary and multicomponent systems. For highly dense systems, for example, typical liquids, it is well known that molecular test insertion methods fail. In this article, we use the Continuous Fractional Component Monte Carlo (CFCMC) method to directly calculate thermodynamic factors by adopting the PWTPI method. The CFCMC method uses fractional molecules whose interactions with their surrounding molecules are modulated by a coupling parameter. Even in highly dense systems, the CFCMC method efficiently handles molecule insertions and removals, overcoming the limitations of the PWTPI method. We show excellent agreement between the results of the PWTPI and CFCMC methods for the calculation of thermodynamic factors in binary systems of Lennard-Jones molecules and ternary systems of Weeks-Chandler-Andersen molecules. The CFCMC method applied to calculate the thermodynamic factors of realistic molecular systems consisting of binary mixtures of carbon dioxide and hydrogen agrees well with the NIST REFPROP database. Our study highlights the effectiveness of the CFCMC method in determining thermodynamic factors for diffusion, even in densely packed systems, using relatively small numbers of molecules.

A framework for subsurface monitoring by integrating reservoir simulation with time-lapse seismic surveys.

Reservoir simulations for subsurface processes play an important role in successful deployment of geoscience applications such as geothermal energy extraction and geo-storage of fluids. These simulations provide time-lapse dynamics of the coupled poromechanical processes within the reservoir and its over-, under-, and side-burden environments. For more reliable operations, it is crucial to connect these reservoir simulation results with the seismic surveys (i.e., observation data). However, despite being crucial, such integration is challenging due to the fact that the reservoir dynamics alters the seismic parameters. In this work, a coupled reservoir simulation and time-lapse seismic methodology is developed for multiphase flow operations in subsurface reservoirs. To this end, a poromechanical simulator is designed for multiphase flow and connected to a forward seismic modeller. This simulator is then used to assess a novel methodology of seismic monitoring by isolating the reservoir signal from the entire reflection response. This methodology is shown to be able to track the development of the fluid front over time, even in the presence of a highly reflective overburden with strong time-lapse variations. These results suggest that the proposed methodology can contribute to a better understanding of fluid flow in the subsurface. Ultimately, this will lead to improved monitoring of reservoirs for underground energy storage or production.

Simulation of the inelastic deformation of porous reservoirs under cyclic loading relevant for underground hydrogen storage.

Subsurface geological formations can be utilized to safely store large-scale (TWh) renewable energy in the form of green gases such as hydrogen. Successful implementation of this technology involves estimating feasible storage sites, including rigorous mechanical safety analyses. Geological formations are often highly heterogeneous and entail complex nonlinear inelastic rock deformation physics when utilized for cyclic energy storage. In this work, we present a novel scalable computational framework to analyse the impact of nonlinear deformation of porous reservoirs under cyclic loading. The proposed methodology includes three different time-dependent nonlinear constitutive models to appropriately describe the behavior of sandstone, shale rock and salt rock. These constitutive models are studied and benchmarked against both numerical and experimental results in the literature. An implicit time-integration scheme is developed to preserve the stability of the simulation. In order to ensure its scalability, the numerical strategy adopts a multiscale finite element formulation, in which coarse scale systems with locally-computed basis functions are constructed and solved. Further, the effect of heterogeneity on the results and estimation of deformation is analyzed. Lastly, the Bergermeer test case-an active Dutch natural gas storage field-is studied to investigate the influence of inelastic deformation on the uplift caused by cyclic injection and production of gas. The present study shows acceptable subsidence predictions in this field-scale test, once the properties of the finite element representative elementary volumes are tuned with the experimental data.

Experimental characterization of [Formula: see text]/water multiphase flow in heterogeneous sandstone rock at the core scale relevant for underground hydrogen storage (UHS).

Geological porous reservoirs provide the volume capacity needed for large scale underground hydrogen storage (UHS). To effectively exploit these reservoirs for UHS, it is crucial to characterize the hydrogen transport properties inside porous rocks. In this work, for the first time in the community, we have performed [Formula: see text]/water multiphase flow experiments at core scale under medical X-ray CT scanner. This has allowed us to directly image the complex transport properties of [Formula: see text] when it is injected or retracted from the porous rock. The important effective functions of capillary pressure and relative permeability are also measured, for both drainage and imbibition. The capillary pressure measurements are combined with MICP data to derive a receding contact angle for the [Formula: see text]/water/sandstone rock system. The rock core sample is a heterogeneous Berea sandstone (17 cm long and 3.8 cm diameter). Our investigation reveals the interplay between gravitational, capillary, and viscous forces. More specifically, it illustrates complex displacement patterns in the rock, including gravity segregation, enhancement of spreading of [Formula: see text] due to capillary barriers, and the formation of fingers/channel during imbibition which lead to significant trapping of hydrogen. These findings shed new light on our fundamental understanding of the transport characteristics of [Formula: see text]/water relevant for UHS.

Geomechanical simulation of energy storage in salt formations.

A promising option for storing large-scale quantities of green gases (e.g., hydrogen) is in subsurface rock salt caverns. The mechanical performance of salt caverns utilized for long-term subsurface energy storage plays a significant role in long-term stability and serviceability. However, rock salt undergoes non-linear creep deformation due to long-term loading caused by subsurface storage. Salt caverns have complex geometries and the geological domain surrounding salt caverns has a vast amount of material heterogeneity. To safely store gases in caverns, a thorough analysis of the geological domain becomes crucial. To date, few studies have attempted to analyze the influence of geometrical and material heterogeneity on the state of stress in salt caverns subjected to long-term loading. In this work, we present a rigorous and systematic modeling study to quantify the impact of heterogeneity on the deformation of salt caverns and quantify the state of stress around the caverns. A 2D finite element simulator was developed to consistently account for the non-linear creep deformation and also to model tertiary creep. The computational scheme was benchmarked with the already existing experimental study. The impact of cyclic loading on the cavern was studied considering maximum and minimum pressure that depends on lithostatic pressure. The influence of geometric heterogeneity such as irregularly-shaped caverns and material heterogeneity, which involves different elastic and creep properties of the different materials in the geological domain, is rigorously studied and quantified. Moreover, multi-cavern simulations are conducted to investigate the influence of a cavern on the adjacent caverns. An elaborate sensitivity analysis of parameters involved with creep and damage constitutive laws is performed to understand the influence of creep and damage on deformation and stress evolution around the salt cavern configurations. The simulator developed in this work is publicly available at https://gitlab.tudelft.nl/ADMIRE_Public/Salt_Cavern .

Pore-scale modelling and sensitivity analyses of hydrogen-brine multiphase flow in geological porous media.

Underground hydrogen storage (UHS) in initially brine-saturated deep porous rocks is a promising large-scale energy storage technology, due to hydrogen's high specific energy capacity and the high volumetric capacity of aquifers. Appropriate selection of a feasible and safe storage site vitally depends on understanding hydrogen transport characteristics in the subsurface. Unfortunately there exist no robust experimental analyses in the literature to properly characterise this complex process. As such, in this work, we present a systematic pore-scale modelling study to quantify the crucial reservoir-scale functions of relative permeability and capillary pressure and their dependencies on fluid and reservoir rock conditions. To conduct a conclusive study, in the absence of sufficient experimental data, a rigorous sensitivity analysis has been performed to quantify the impacts of uncertain fluid and rock properties on these upscaled functions. The parameters are varied around a base-case, which is obtained through matching to the existing experimental study. Moreover, cyclic hysteretic multiphase flow is also studied, which is a relevant aspect for cyclic hydrogen-brine energy storage projects. The present study applies pore-scale analysis to predict the flow of hydrogen in storage formations, and to quantify the sensitivity to the micro-scale characteristics of contact angle (i.e., wettability) and porous rock structure.